How to look for supersymmetry under the lamppost at the LHC 
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We apply a model-independent, agnostic approach to the collider phenomenology of supersymme- 
try (SUSY), in which all mass parameters are taken as free inputs at the weak scale. We consider 
the gauginos, higgsinos, and the first two generations of sleptons and squarks, and analyze all pos- 
sible mass hierarchies among them (4x8! = 161, 280 in total) in which the lightest superpartner is 
neutral, leading to missing energy. In each case, we identify the full set of the dominant (i.e. least 
suppressed by phase space, small mixing angles or Yukawas) decay chains originating from the 
lightest colored superpartner. Our exhaustive search reveals several quite dramatic yet unexplored 
multilepton signatures with up to 8 isolated leptons (plus possibly up to 2 massive gauge or Higgs 
bosons) in the final state. Such events are spectacular, background-free for all practical purposes, 
and may lead to a discovery in the very early stage (~ 10 pb _1 ) of LHC operations at 7TeV. 

PACS numbers: 14.80.Ly,12.60.Jv,13.85.-t 



The ramping operations at the CERN Large Hadron 
Collider (LHC) have begun the long awaited and historic 
exploration of the TeV scale, where new physics beyond 
the standard model (SM) is expected to emerge. Among 
the multitude of new scenarios, low energy supersym- 
metry (SUSY) has long been the primary target of the 
LHC experiments, not just because it is well motivated 
theoretically but also because its generic discovery 
signatures cover a much wider class of models 0] . 

By itself, SUSY is very predictive, as it fixes the spins 
and couplings of the new particles (the superpartners) . 
Unfortunately, this is not sufficient to pin down its pre- 
cise collider discovery signatures, as the latter crucially 
depend on the SUSY mass spectrum, which is in turn 
determined by the mechanism of supersymmetry break- 
ing. Alas, almost 40 years of model building effort since 
the discovery of supersymmetry have failed to produce a 
single, universally accepted model of SUSY breaking. 

Given one's utter ignorance about the expected pat- 
tern of SUSY masses, in this letter we adopt a most con- 
servative, agnostic approach, where the masses of all su- 
perpartners are treated as free inputs at the weak scale. 
We shall then consider all possible hierarchical patterns 
among them, and identify the set of dominant (in the 
sense defined below) collider signatures in each case. In 
our quest for interesting models, we shall be guided by ex- 
perimental pragmatism instead of theoretical prejudice. 

The purpose of this study is twofold. First, most pre- 
vious collider studies of SUSY have been performed for 
specific SUSY benchmark points, typically chosen within 
some minimal model such as "minimal supcrgravity" 
(MSUGRA) 0. We will therefore be interested in un- 
covering new types of signatures which may have been 
missed in the standard benchmark approach. Secondly, 
we shall focus our search on signatures with a high num- 
ber of isolated leptons, which constitute the proverbial 
"smoking gun" for new physics. For example, the inclu- 
sive trilepton channel is already recognized as "the golden 
mode" for an early SUSY discovery at hadron colliders. 
One of our main results here will be the identification of 



TABLE I: The set of SUSY particles considered in this anal- 
ysis, shorthand notation for each multiplet, and the corre- 
sponding soft SUSY breaking mass parameter. 
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a number of new SUSY mass patterns whose dominant 
signatures have up to eight leptons in the final state. 

Our setup is as follows. We take the usual superpartner 
content of the Minimal Supcrsymmetric Standard Model 
(MSSM) listed in Table Q] For simplicity, in this letter 
we shall consider just two degenerate light generations 
of sfermions. Third generation effects can be trivially 
incorporated in the discussion Q, and only complicate 
the bookkeeping. Given the 9 input mass parameters in 
Table HI in general there are 9! = 362, 880 possible or- 
derings among them, each leading to a distinct pattern 
(hierarchy) of sparticlc masses. We shall use the short- 
hand notation from Table Q] to label each hierarchy: for 
example, GQUDHLWEB is a model with M G > M Q > 
Mu > M D > M H > M L > M w > M E > M B . 

Our first goal will be to identify the main collider sig- 
natures for each hierarchy. As in any discussion on SUSY 
collider phenomenology, our starting point is the fate and 
then the nature of the Lightest Supersymmetric Particle 
(LSP), which we shall generically denote by C. For our 
main analysis, we shall assume that i?-parity is conserved 
(or very weakly broken), so that C is stable on the scale 
of the detector. (We briefly discuss the i?-parity violating 
option at the end.) Then, the original 9! model hierar- 
chies can be classified into the following three categories: 

I. CHAMPs. In the 8! = 40, 320 cases with C = E, 
the LSP is an electrically charged, color-neutral parti- 
cle (the right-handed slepton Br). The corresponding 
generic collider signature is a long-lived charged massive 
particle (CHAMP) [f|, regardless of the particular order- 
ing of the heavier sparticles. 
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FIG. 1: Graphical representation of the allowed transitions 
between the SUSY states from Table [I] One (two, three) 
parallel lines represent two- (three-, four-) body decays. The 
relative suppression of each decay mode is indicated by the 
line type. The identity of the resulting SM decay products is 
denoted by the line color: red for a jet j, blue for a lepton I 
and green for a massive boson v = {W ± , Z, h} (which may 
be either on-shell or off-shell). 

II. R-hadrons. In 4 x 8! = 161,280 of the remaining 
hierarchies C G {Q, U, D, G}, the LSP is a colored parti- 
cle, and the generic searches for stable R-hadrons apply 

Again the ordering of the heavier particles is not 
particularly important. 

III. Missing transverse energy. In the remaining 4 x 
8! = 161,280 cases C G {L,B,W,H} and the LSP is a 
weakly-interacting, electrically neutral particle. Its pro- 
duction will lead to missing transverse energy ($t) in the 
detector. Now, however, the signatures crucially depend 
on the ordering of the heavier particles, since it is not 
feasible to look for $t inclusively. Our goal here is to 
fully classify these 161, 280 models according to their col- 
lider phenomenology. Unlike previous general approaches 
0] , which employed scans of the multi-dimensional SUSY 
parameter space, here we would like to avoid scanning, 
keeping the discussion simple and qualitative. 

Both of the currently operating high energy colliders 
(the Tevatron at Fermilab and the LHC at CERN) are 
hadron machines, at which the total production is ex- 
pected to be dominated by the strong production of col- 
ored supcrpartners. Correspondingly, the starting point 
for our classification will be the nature of the lightest col- 
ored supcrpartner (LCP), denoted by C. Then, each of 
the 161,280 missing energy hierarchies at hand can be 
represented by a particular ordering 

x...xCy...yC, (1) 

where the it's stand for inconsequential entries, C G 



TABLE II: Number of hierarchies for the various dominant 
decay modes of the LCP C. 
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{G, Q, U, D}, y G {L, B, W, H, E} and C G {L, B, W, H}. 
The dominant collider signature for each model hierarchy 
([Tj) will be determined by the inclusive pair production 
of C and its dominant subsequent decays. 

Our key idea here is that once a given hierarchy ([TJ is 
assumed, the dominant decay modes of C are uniquely de- 
termined, since supersymmetry predicts all supcrpartner 
couplings. In our analysis, we shall assume that there 
are no accidental phase space suppressions due to any 
two mass parameters from Table [I] being very close. This 
assumption also guarantees that the chargino and neu- 
tralino mixing angles are small and the mass eigenstatcs 
arc roughly aligned with the interaction eigenstatcs. One 
can then use the simple chart in Fig. [T] to identify the 
dominant (i.e. least suppressed) decay modes of C, which 
we label by the number of leptons n£ (blue lines), number 
of jets nj (red lines) and number of massive bosons n v 
(green lines) encountered along the way. Solid lines in the 
figure correspond to 2-body decays which do not suffer 
from any (chargino or ncutralino) mixing angle suppres- 
sion (MAS); dashed lines indicate either 2-body decays 
with MAS or 3-body decays with no MAS; and finally, 
dotted lines stand for either 3-body decays with MAS or 
4-body decays with no MAS. We then count the num- 
ber of mass hierarchies (UJ) which exhibit a dominant de- 
cay channel for C with a given set of (ri£, n„, rij), and 
show the result in Table QTJ At times, there can be sev- 
eral dominant decay modes of C. For example, consider 
xxxxQW BLH . One can get (ni,n v ,nj) = (0,1,1) in 
two ways: Q — s- W — > H or Q —> B — > H. It is also pos- 
sible to have (rii, n v , rij) = (2, 0, 1) in two different ways: 
Q^W^L^H or Q^B^L^H. Therefore, 
xxxxQW BLH contributes one entry to each of the two 
boxes (ri£,n v ,rij) = (0,1,1) and (ne,n v ,rij) = (2,0,1) 
in Table |ITJ As a result, the total number of entries 
(203,184) in Tabic [TTI is larger than the total number of 
hierarchies (161,280). 

Table HT1 leads to some interesting conclusions. For ex- 
ample, we see that the purely hadronic signatures of rij 
jets and fr alone cover a very large fraction (~ 65%) of 
all possible SUSY hierarchies with a neutral LSP. There 
is also a sizable fraction of models which can be explored 
via the standard searches for signatures with one, two 
or three leptons. Keep in mind that the LCP's are pro- 
duced in pairs, so the collider signature is obtained by 
doubling the number of leptons and jets displayed in the 
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TABLE III: Number of hierarchies for the maximally leptonic 
decay modes of the LCP C. 
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table and is given by 2ng + 2n v + 2rij + fr- Perhaps 
the most important result from Table [TTI is that there is a 
non-negligible fraction of SUSY models (~ 1%) in which 
one of the dominant decay modes of the LCP yields 4 iso- 
lated leptons, and the corresponding collider signature is 
8 leptons plus 

We now repeat the same analysis, only this time from 
the set of all dominant decay modes of C, we select the 
one with the largest m, and in case of a tie for m, we 
pick the chain with the larger n v . We refer to such decay 
modes of the LCP as "maximally leptonic". The new 
tally is displayed in Table IIII1 where now the sum of all 
entries equals the total number of signatures 161,280. 

It is instructive to apply our general formalism to the 
familiar MSUGRA model, where the soft SUSY mass 
parameters of Table U have common values (Mo for the 
scalar supcrpartners and M1/2 for the gauginos) at the 
grand unification scale. First, one may ask how many of 
the 161,280 missing energy hierarchies are actually rep- 
resented in MSUGRA. The answer is provided in Fig. [5J 
in which we divide the usual (Mo, Mi 12) plane into dis- 
joint color-coded areas, according to the observed mass 
pattern (fTJ) . The allowed region contains only 47 diffcr- 
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FIG. 2: A slice through the MSUGRA parameter space for 
fixed Ao = 0, tan ft = 10 and yi > 0. Yellow shaded regions 
are ruled out by direct searches at LEP or by requiring a neu- 
tral LSP. The remaining area is color-coded according to the 
sparticle hierarchy type JT]). The solid black lines delineate 
regions with the same maximally leptonic C decay mode. 



TABLE IV: Input soft SUSY mass parameters (in GeV) for 
the xxGQWLBEH study points used for Figs. [3] and 0] 
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ent hierarchies (some areas are too small to be readily 
noticeable with the naked eye). Comparing this to the 
total number of 161,280 possibilities, one gets an idea of 
the limitations of MSUGRA as a benchmark scenario. 

We also use Fig. [2] to illustrate the extent to which 
the MSUGRA model is able to cover the generic SUSY 
signature space. The black solid lines in Fig. [2] divide 
the allowed portion of the (M ,M 1 / 2 ) plane into dis- 
joint regions, classified according to the maximally lep- 
tonic C decay mode, labelled by (ni,n v ,rij). We find 
that MSUGRA exhibits only 4 out of the 26 possibil- 
ities found in Table IIIII At small values of Mo, the 
LCP is a right-handed squark (either Hr or dp), whose 
single dominant decay mode is directly to the bino-like 

LSP: {uR,dn} -4 b°. As the value of Mq increases 
and the squarks get heavier, the gluino eventually be- 
comes the LCP and its maximally leptonic decay mode 

is g ^4 w° -A b°. Upon further increasing Mo, we even- 
tually enter the focus point region ||, where at first 



Mw > Mh > Mb and the winos decay as { 



,r,± ,r,o 



} A 



{h^, h^, h d } A b°, while in the rightmost portion of the 
plot one finds Mw > Mb > Mh with the direct transi- 
tions {w ± ,w } A {h ± , h°, hP d }. Of course, Fig. [2] only 
displays dominant LCP decay modes. Some parame- 
ter points in MSUGRA may contain other, longer decay 
chains with a higher number of leptons, but those would 
be subdominant and therefore suppressed by branching 
fractions. For example, the hierarchy xxxU HLWEB has 

a 4-lepton chain ur ^> h° —t et, — > w° —t b°, but it 

has two bottlenecks at ur A h° and h° —> e L , see Fig.Q] 
In the rest of this letter we study in detail one 



example 11| of a maximally leptonic decay chain with 4 
leptons, where the corresponding collider signal is 8 iso- 
lated leptons plus jets plus missing energy. For concrete- 
ness, consider the hierarchy xxGQWLBEH at several 
different study points, defined in Table IIVI and chosen 
"under the lamppost" , i.e. to maximize the 8-lepton sig- 
nal rate. For this hierarchy, 8 lepton events arise from the 
inclusive pair production of left-handed squarks ul^l, 
followed by 



ul,u l -> w 



II A6 u ^Ky d . (2) 

Ignoring phase space suppression factors, third genera- 



4 




300 400 500 600 700 800 900 

M„(GeV) 

FIG. 3: Inclusive cross-sections (in fb) at the LHC for 8 lepton 
events from the exclusive chain in for the study points from 
Table HV1 (solid lines); inclusive cross-sections for left-handed 
squark pair production (dotted); and cumulative branching 
fraction of the 8 lepton final state (dashed, read on the right). 



tion effects and the masses of the W^, Z and ft,, it is easy 
to estimate the individual branching ratios for this chain 
as follows 3: BR(g ->■ q L + j) = 1; BR{q L w° + j) ~ 
|; BR(w° -> £± + £T) ~ |; BR(£jr -> b° + £±) = 1; 
BR(b° I| + fF) ~ |; ft° /d + £±) = 1. 

Multiplying these results and squaring, we expect the 
cumulative branching fraction for 8 lepton events to be 
around 0.7%. This is confirmed by the dashed line in 
Fig. [31 which approaches this asymptotic value for large 
mass splittings (large Mq). We see that to maximize the 
overall branching fraction, we must scale the SUSY mass 
spectrum up. Yet to maximize the production rate, we 
must scale the spectrum down. An optimum compromise 
is therefore found for intermediate values of the SUSY 
masses, as shown in Fig. [3l The study points in Table |IV] 
were picked by varying Mq, fixing Mq = Mq + 100 GeV, 
and choosing the rest of the spectrum from a coarse scan 
to maximize the 8 lepton rate shown by the solid lines 
in Fig. [3] These study points were then processed with 
the PYTHIA @ event generator and the PGS [T3] de- 
tector simulator for the case of an LHC at 7 TeV center- 
of-mass energy and just 1 fb -1 of data. Because these 
events are essentially background free, we simply count 



PGS-reconstructed isolated leptons with default pr cuts 
of 3 GeV for muons and 10 GeV for electrons, and display 
the result in Fig. 01 Due to the imperfect detector accep- 
tance, we got only one 8 lepton event in 1 fb -1 . Never- 
theless, for Mq < 500 GeV typically there are a handful 
of 7 lepton events and hundreds of 4 lepton events, which 
are already very clean. Since only 5 clean events are suf- 
ficient for discovery, the region Mq < 500 GeV can be 
probed with as little as 10 pb _1 of data. 

We note in passing that one may expect even more lep- 
tons if i?-parity is violated and the LSP decays promptly 
to SM particles. For example, lepton number violation of 
the type LLE would lead to 2+2=4 additional leptons, 
bringing the maximum lepton total per event to 12. 

In conclusion, we have shown that already with its first 
10 pb _1 of data the LHC can start probing SUSY, pro- 
vided one looks in the right place. To this end, it is 
imperative to search for the model hierarchies that are 
most likely to be discovered first. We believe this letter 
outlines the most general strategy to do just that. 
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FIG. 4: (Stacked) number of multi-lepton events in 1 fb -1 
of data at a 7 TeV LHC, for the xxGQWLBEH sparticle 
hierarchy with the mass spectra shown in Table ITVl 
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